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We study a new type of Fermi polaron induced by an impurity interacting with an ultracold Fermi 
superfluid. Due to the three-component nature of the system, the polaron can become trimer-like 
with a non-universal energy spectrum. We identify multiple avoided crossings between impurity- 
and trimer-like solutions in both the attractive and the repulsive polaron spectra. In particular, the 
widths of avoided crossings gradually increase as the Fermi superfluid undergoes a crossover from the 
BCS side towards the BEC side, which suggests instabilities towards three-body losses. Such losses 
can be reduced for interaction potentials with small effective ranges. We also demonstrate, using the 
second-order perturbation theory, that the mean-field evaluation of the fermion-impurity interaction 
energy is inadequate even for small fermion-impurity scattering lengths, due to the essential effects 
of Fermi superfluid and short-range physics in such a system. Our results are practically useful for 
cold atom experiments on mixtures. 


I. INTRODUCTION 

Quasiparticles serve as the cornerstone of complex col¬ 
lective phenomena in interacting many-body systems. As 
an example, the polaron, an impurity dressed by particle- 
hole excitations from the environment, is a typical quasi¬ 
particle that has stimulated much interest in the study of 
solid-state and cold-atom systems with highly polarized 
components. In cold atoms, polarons have been success¬ 
fully explored in both the attractive and the repulsive 
branches of impurities interacting with a Fermi sea of 
identical fermions m- These studies are crucial for un¬ 
derstanding the nature of normal states either against 
pairing physics in the case of strong attraction [5lH0], or 
in the context of itinerant ferromagnetism with strong re¬ 
pulsion mm- However, much less is known about the 
fate of impurities immersed in a Fermi superfluid, where 
the impurity can be dressed by pairs of particles in the su¬ 
perfluid M- Given the growing number of experiments 
on atomic mixtures, it becomes pressing to understand 
the fundamental physics of the underlying quasiparticles, 
such as polarons, in these systems. 

Particularly, in the recent ENS experiment, a mixture 
of Bose and Fermi superfluids has been realized [15]. The 
two superfluids not only coexist, but also interact with 
each other, constituting the most charming character of 
this system. How to characterize the interaction energy 
between these two superfluids, however, is still an open 
question. Exact three-body calculations show that, in the 
BEC limit of fermions, the interaction between a bosonic 
atom and a molecule of two fermions cannot be faith- 
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fully described by the mean-field theory, even when the 
boson-fermion scattering length is small mm- It then 
becomes essential to investigate the validity of the mean- 
field theory on the many-body level throughout the whole 
BCS-BEC crossover regime of fermions, and identify the 
underlying mechanism wherever it fails to apply. Here, 
we approach this goal via the study of polarons. 

We use the variational approach to investigate the po¬ 
laron physics when an impurity is immersed in a Fermi 
superfluid with tunable interactions. With a pairing su¬ 
perfluid as the background, the polaron wave function 
naturally acquires trimer-like terms which originate from 
the impurity-induced pair-breaking processes in the su¬ 
perfluid bath. As a result, the polaron spectrum is 
composed of both impurity- and trimer-dominated so¬ 
lutions. Due to the couplings inbetween, many avoided 
level crossings occur, whose widths depend sensitively on 
both the interaction strength and the effective range of 
the system. In particular, as fermions are tuned from the 
BCS to the BEC side, the widths of the avoided crossings 
gradually increase, which suggests instabilities towards 
three-body losses. On the contrary, the widths become 
narrower at smaller effective ranges, which can help sta¬ 
bilizing the system. For small fermion-impurity scatter¬ 
ing lengths, the second-order perturbation theory cap¬ 
tures the essential effects of Fermi superfluid and short- 
range physics in the system, and consequently shows the 
insufficiency of mean-field descriptions. These results 
rectify our understanding of interaction effects in multi- 
component mixtures on the mean-field level, and provide 
guidance for maintaining stability from three-body losses 
in cold atom experiments on mixtures. 

The paper is organized as the following: in Sec. [H| 
we outline the variational approach adopted in this work 
and provide details on the derivation of the equations 
used for the numerical evaluation of the polaron energy. 
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In Sec. m we present the full polaron spectra and the 
polaron residue. We then apply the second-order pertur¬ 
bation theory in Sec. 23 and discuss the insufficiency of 
the Hartree-type mean-field characterization of the po¬ 
laron energy. In Sec. [V] and Sec. we study the effects 
of fermion-fermion interaction and the effective range on 
the polaron spectrum, respectively. We examine the po¬ 
tential polaron to molecule transition in the system in 
Sec. |VII| Finally, we summarize in Sec. VIII 


Here, is the annihilation op¬ 

erator for the Bogoliubov quasiparticles, where 77 ^ = 
— 77 ^ = 1 , 77 k = yj (FV + e k — aO/ ( 2 -Ek), and ^k = 
— 77 k ; the corresponding quasiparticle energy E\^ = 
yf (ek — /i ) 2 + A 2 with A the pairing order parameter. 

To depict the impurity-induced polaron excitations in 
a Fermi superfluid, we adopt an ansatz inspired by that 
used in a highly-polarized two-component Fermi gas [5j 

EDIEU: 


II. MODEL 

We start from the Hamiltonian of our system JC = 
H — fiNp {Np is the total number of fermions): 

IC = ka + y a kt a q-k4. a q- k '-l- a k't 

k,cr k,k',q 

+ ^V k V k +^7 V V aigbq-kbq-vava, (1) 

k k,k',q cr 

where a\^ a and 6k are respectively the annihilation oper¬ 
ators for the superfluid fermions and the impurity atom 
with the dispersion ek = k 2 /(27n) (h is taken to be 
unity); g is the chemical potential of the two-species 
(cr =t,|) spin-balanced Fermi gas. gft (g&) is the bare 
fermion-fermion (fermion-impurity) interaction, which is 
related to the scattering length (clr) via the stan¬ 
dard renormalization relation: PS l /gp = m/ (47t ap) — 
1/V l/(2ck), where V is the quantization volume and 

f3 = ff, fi. For simplicity, we only consider the case where 
the impurity has the same mass as that of a fermion, 
and interact equally with two fermion species. Our re¬ 
sults can be straightforwardly generalized to cases with 
unequal masses or imbalanced interactions. 

The Fermi superfluid at zero temperature can be de¬ 
scribed by the standard BCS wave function m : 

|BCS) = H(«k + ^ka k t a - k 4.)l vac } ~ qa_ kJ .a kt |vac). 

k k 

( 2 ) 


I p )q = ( V’q^q + V^k'^Q-k'-kd'i-di) l BCS ), 

V kk' / 

(3) 

where Q indicates the center-of-mass momentum of the 
polaron. The second term in the bracket, which effec¬ 
tively describes impurity-induced pair breaking in the 
superfluid and is therefore trimer-like, includes contri¬ 
butions from excitations like a k/ ^a_k^, aj^a-kq? a k't a H 
or a_k^U-kt- As a first attempt at the problem, we only 
keep excitations to the lowest order [22]. The existence of 
trimer-like terms in polaron wave functions is unique for 
a pairing-superfluid background, which significantly af¬ 
fects the polaron spectra. Trimers in a many-body back¬ 
ground have received much attention recently in different 
contexts [231427] , and ours is a new platform where such 
few-body effects can be observed. 

The ground state solution can be obtained by minimiz¬ 
ing E p = (P\JC\P) - Ebcs , with E B cs = 2k( e k — g — 
.Ek) + Xlk A 2 /(2Ek). Note that it is convenient to ex¬ 
press the fermion-fermion interaction term in 1C in terms 
of the scattering between different quasiparticle states: 

k, q — k —>> k 7 , q — kb Since in the large- k limit, 

77k —> 1, 77k 0, only one term has finite contribution 

to the energy: 

9ft t t (A\ 

y / J ^k^q—k^q—k'^k'C^kt^q—k^^U - • ( 4 ) 

q,k,k' 


We then have the equations: 


{Ep - e Q )^Q = y ( 2^ \v^\ 2 l/jQ + ^2 ^k'V’kk' + ^^k'^kk' 


(5) 


kk' 


kk' 


Akk'^kk' = yr f ^k' ^ ^k"V’kk" “ ^k ^ + ^k ^ Tik'^k^k' “ ^k' ^ ^k'^kk" + ^k^k'^Q + ^k^k'^Q 

V k" k" k" k" j 


9ft / 

— 77k^k' / J Vk // ,k+k / —k 


/' 77k" ^k+k' —k" • 


(6) 


As gp {/3 = ff, fi) would become vanishingly small after renormalization, terms including gp^2 k ^k --- 5 where 77k 


rsj 
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1/k 2 at large momentum fc, should also vanish. We therefore neglect these terms and define: 

Ak = fffi^k^Q + y> k ^,kQ; (7) 

k 7 

^k = fffit^k^Q + y^^^,k); (8) 

k 7 

Ck = gft ^^k'^-k'^k'^k-k'* (9) 

k 7 


From these, we can get V’Q and V^k 7 in te rms of 
Ak,£?k,Ck, which, when plugged into Eqs. <om>> yield 
a set of closed equations: 


w ^kk 7 
k 7 


Ay, = 


(2v k v k 'Av u k u k >A k > 


-E 


E ^ Ep cq 

^ ^fe^fe/C'k+k 7 


A 


kk 7 


^fe^k-k 7 

Ak 7 ,k-k 7 


^kk 7 

ri sr~^ ^k—k 7 Uk' Ay,/ 

Ck = E ^ k , k _ k , 

h ' ’ 


( 10 ) 


( 11 ) 


with Ak,k 7 = E p — 6q— k—k 7 — Ek — Ek'. Physically, as 
we only consider the case of spin-independent impurity- 
fermion interactions, we expect Ak = By,, ^k,k 7 = Vv,k- 
We have numerically checked that this is the case. 

It is also important to note that Eqs. Pint are able 
to reproduce the exact three-body equations when the 
fermion density is sent to zero. While with a finite 
fermion density, the interplay of Fermi superfluid and 
trimer physics gives rise to intriguing polaron properties 
as discussed below. In practice, E p and variables Ak, Ck 
can be numerically solved by imposing a certain momen¬ 
tum cutoff k c . Physically, k c corresponds to setting a 
finite effective range ro ^ 1 Jk c in the two-body collision. 
The coefficients in ansatz (|3| can be easily deduced from 
the Schrodinger’s equation §, which gives the polaron 
residue [28]: 


Z = 






Q 


+ X)k,k' ^kk' 


( 12 ) 


Apparently, the polaron residue Z r\j 1 (Z rsj 0) when the 
system is impurity (trimer) dominated. In the rest of the 
work, we will focus on polarons in the Q = 0 sector. 



FIG. 1: (Color Online) Polaron spectrum (a) and residue 
(bl,b2) as functions of the fermion-impurity interaction 
strength. The fermion-fermion interaction is at resonance 
(as = oo). (a) Blue solid, red dashed and black solid lines 
correspond to the lowest, second lowest and higher branches 
of the spectrum. Red dots and magenta dash-dotted line 
show, respectively, the perturbative energies Ept/Ef and 
(E pt + Ep T )/Ef. (bl,b2) The polaron residue for the two 
lowest branches (bl) and several higher branches (b2). The 
magenta dash-dotted curves in (bl,b2) are the residue results 
based on the perturbative wave function. The cutoff momen¬ 
tum k c = 10/cf, and the unit of energy Ef — k 2 F /2m. 


III. POLARON ENERGY AND RESIDUE 


By explicitly solving Eqs. dioITTI ), we get the typi¬ 
cal polaron spectrum for a unitary Fermi superfluid (see 
Fig. [lja)). The spectrum exhibits many avoided level 
crossings between impurity- and trimer-dominated so¬ 
lutions, which are also apparent in the residue plots in 
Fig. [ljbl,b2). Particularly, a wide avoided crossing ex¬ 
ists between the two lowest branches at (&F&fi) _1 ^ —3, 


where kp is the Fermi momentum of a non-interacting 
Fermi gas with the same number density as the Fermi 
superfluid. When a^ 1 is tuned through the avoided cross¬ 
ing, the wave function continuously evolves from impu¬ 
rity dominated (Z ~ 1) to trimer dominated (Z ~ 0) 
in the lowest branch [29]. Similar atom-trimer crossover 
has been reported for the ground state of a two-channel 
model in Ref. |14], where equal fermion-fermion and 
fermion-impurity interactions are considered. A polaron 
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(a) 



(b) 


due to the enhanced coupling between the impurity and 
the trimer terms in Eq. ©• This repulsive branch even¬ 
tually runs into a dense spectrum of trimer-dominated 
solutions. 




^ 


+ ... 


FIG. 2: (Color Online) Diagrams for the second-order per¬ 
turbative corrections in |afi|: (a) and (b) are, respectively, 
without and with the intermediate scatterings between su¬ 
perfluid fermions. The black solid (red dashed) line is the 
propagator for impurity (Bogoliubov quasiparticles). The 
hollow (solid) red circles represent the bare (renormalized) 
fermion-fermion interactions, and the black circles represent 
the fermion-impurity interactions. 


to molecule transition prroim] can also occur within the 
second lowest branch of our system, when the fermion- 
impurity interaction is tuned towards the BEC limit. 

Compared to the two lowest branches, the higher 
branches at positive a g in Fig. [lja) show more interest¬ 
ing properties, as multiple avoided crossings occur with 
widths tunable by a&. At small ag, the avoided cross¬ 
ings are very narrow, which allows us to identify a re¬ 
pulsive atomic branch by following the trajectory of the 
impurity-dominated solutions (Z ~ 1) while a a varies. 
The energy of the repulsive branch increases with ag, 
while the width of avoided crossings becomes broader 


IV. PERTURBATIVE CORRECTIONS 


To gain further insights into the polaron state, we ap¬ 
ply the second-order perturbation theory at small |afi|. 
Intuitively, up to a|, the perturbative energy caused by 
the fermion-impurity scattering can be written as: 


Epx = 


4tt a fi p +2 /^fij2V 
m mV ' ^ 

kq 


2 ( 1 

M2 Vk 




6— q_ k T Eq T Efc 


■ E n + Eh 


(13) 


Here, the first term is the mean-field contribution, with 
p the total density of fermions; while the second term 
accounts for the intermediate scattering processes shown 
in Fig. |2^a). Eq. (13) can be continuously reduced to 


the case of non-interacting fermions m with decreas¬ 
ing fermion-fermion interaction. However, for a uni¬ 
tary Fermi superfluid, Fpx is well above the polaron 
energy E p even at relatively small |afi| (see Fig. Oja)). 
This can be attributed to the effects of intermediate 
fermion-fermion scatterings, as shown diagrammatically 
in Fig. |2|b). The sum of all relevant diagrams gives an 
additional energy correction up to the order of a|: 




(V'k'Vq T r UqVk')T\ i ,q;k' ,k+q—k' {'Uk'V k+q—k' t%%+q-k') 


mV 


k,q,k' 


(C —q_k T -F/c T Eq){€ — q—k T Efc 


■E } 


^k^k',k+q—k' — ^k^q^k'^ k+q—k' 


mV 

47raff 


-E^ + E 


k+q-k'J 


U k" U k+q-k" 


(14) 


k 26k q—k T E]^n T F/k+q— k" 


As shown in Fig. [lja), the inclusion of Ep T consider¬ 
ably lowers the perturbative energy based on Eq. (13). 
The remaining discrepancy is attributed to higher-order 
impurity-fermion scattering processes beyond the dia¬ 
grams in Fig. [2j which lead to energy corrections on the 
order of with n > 3. 

The second-order perturbation results reveal several 
remarkable properties of polarons in the presence of pair¬ 
ing superfluid. First, through the inclusion of Fp T (Eq. 
©) and u k , v k , E k (Eqs. ( Pim p , the second-order 
terms rely crucially on the interaction between fermions. 
Thus, sweeping the Fermi superfluid across the reso¬ 
nance is expected to significantly affect the polaron. Sec¬ 
ond, the summations in Fpx and Ep T both scale loga¬ 
rithmically with the momentum-cutoff k c (see Fig. [ 3 ]). 


This scaling relation indicates that E p is generally non- 
universal, and that the mean-field evaluation of E p is 
inadequate even for small |afi|. The reason of the non- 
universality can be attributed to the three-component na¬ 
ture of the system, where the high-energy (short-range) 
detail, or the effective range, plays an essential role. In 
the following, we will show the dramatic effects of the 
fermion-fermion interaction and the effective range on 
the polarons. 


Finally, the polaron residue Z can also be estimated 
based on first-order (in \a^\) wave functions derived 
from the perturbation theory by setting Tq « 1 (see 
Fig.0bl,b2)). 
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FIG. 3: (Color Online) The cutoff-momentum dependence 
of the second-order perturbative corrections: (a) is for Ep,p 
(second-order terms in Eq. (13)); (b) is for Ep T . Here we 
choose interaction parameters (^FUff) -1 = 0.1, (/cf^r) - 1 = 2. 
The unit of energy is taken to be Ef . 


FIG. 5: (Color Online) Effect of effective range on the po- 
laron spectra for the lowest two branches (a) and the higher 
branches (b). The fermion-fermion interaction is at resonance 
(as = oo). The solid and dashed lines are respectively for 
k c = 10 kF and 20 /cf- 



FIG. 4: (Color Online) Effect of fermion-fermion interac¬ 
tion on the polaron energy ( E p ) for the lowest two branches 
(a,b) and higher branches (c,d). (a) E p as functions of 

{kFd a) - 1 , with (/cFUff) -1 = 0 (solid), 0.5 (dashed), and 
— 1 (dash-dotted), (b) E p as functions of (/cf^af) - 1 at a 
fixed (/cFUfi) -1 = —3. (c) E p as functions of (kFd a) - 1 with 
(/cFUff) -1 = 0 (solid), 0.5 (dashed), and —0.3 (dash-dotted), 
(d) E p as functions of (kFd ff) -1 at a fixed (kFd a) - 1 = 
2. In (b,d), the red dots and the magenta dashed-dotted 
line are respectively the perturbative energies Ept/Ef and 
(E pt + E'p T )/Ef . The cutoff momentum k c = 10&f- 


which we attribute to the enlarged phase space of the 
fermion-impurity scattering: as the fermions are more 
tightly paired, the Fermi surface becomes more smeared 
out, and the scattering phase space is less affected by 
the Pauli principle. Accordingly, more impurity-induced 
excitations emerge in the Fermi superfluid, which effec¬ 
tively enhances the impurity-trimer coupling in Eq. <§• 

In addition, from Fig.[4jb,d) one can see increasing de¬ 
viations from the perturbative results when the fermions 
are tuned from the BCS to the BEC side. This suggests 
that a strong Fermi superfluid can dramatically affect 
the polaron spectra through higher-order scattering pro¬ 
cesses, which are beyond the ones shown in Fig. [2] This is 
consistent with the exact calculation of atom-dimer scat¬ 
tering length in the deep BEC limit of fermions maun, 
which includes all orders of scattering processes and 
shows the breakdown of the mean-field prediction. 


VI. EFFECT OF EFFECTIVE RANGE 


V. EFFECT OF FERMION-FERMION 
INTERACTION 


We study the variation of the polaron spectrum as 
the fermion-fermion interaction changes. The results 
are shown in Fig. [4j For both the lowest two branches 
(Fig. B a,b)) and the higher branches (Fig. [4jc,d)), the 
avoided crossings become broader as fermions are tuned 
towards the BEC side. As a result, the repulsive atomic 
branch becomes difficult to identify even for small pos¬ 
itive an (Fig. Qd)), which suggests instabilities towards 
three-body losses. 

A broad avoided crossing implies an enhanced coupling 
between the impurity and the trimer terms in Eq. 


As the effective range (r o) corresponds to the inverse 
of the cutoff momentum (k c ) in our formalism, its effect 
can be studie d by changing k c in numerical simulations 
of Eqs.( Ib|ll ). In Fig. i we show the polaron spectra for 
two different k Cl where Fig.[5ja) and (b) are respectively 
for the lower and the higher branches. A common feature 
is that for larger k C: or smaller ro, the avoided crossings 
move towards weaker fermion-impurity interactions, i.e., 
towards smaller \an\. Accordingly, the coupling between 
impurity- and trimer-dominated solutions also becomes 
smaller, which naturally leads to narrower avoided cross¬ 
ings. In practice, these results suggest that a system with 
a smaller ro can be more stable, as a narrower crossing 
with trimer-dominated solutions makes the decay into 
deep trimers less likely. 
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VII. MOLECULAR STATE AND 
POLARON-MOLECUE TRANSITION 


The general variational wave function for a two-body 
state involving the impurity atom can be written as: 

\ m )q = E (WEkdt + l BCS >- (! 5 ) 

k 

In the case of a spin-independent fermion-impurity inter¬ 
action, the bound state wave functions and are 
decoupled. There is then simply a two-fold degeneracy 
in the two-body sector. We therefore write the ansatz as: 

\M) Q =J2^ Q -A\BCS). (16) 

k 

From the Schrodinger’s equation, the wave functions 
(£>kt should satisfy 

Ecp kf = e Q-k^kt + E k if\^ 

+ y I ttk E M k'</^k't - v k E V k'</?k' + 2 E l w k' I Vk ] , 

V k' k' k' / 

(17) 


where the last term is the Hartree term. As the Hartree 
terms here only serve to shift the molecular energy by a 
vanishingly small amount after renormalization, they are 
not important for the two-body bound states. Similarly, 
we may drop the term g^v k 'Gc'^k' • This leads to the 

equation in the two-body bound state sector: 


V _ _ |^k | 2 _ 

Qfi k Em e Q —k Tk 


(18) 


With the zero energy reference chosen to be the BCS 
ground state energy, the two-body bound state threshold 
F?th in this sector is the quasi-particle excitation gap: 
F?th = A when fi > 0 and E t h = \//S 2 -\- g 2 when 
l± < 0. In Fig. [6|a), we show the lowest two branches 
of the polaron energy and the molecular energy rela¬ 
tive to the threshold. While there are no polaron to 
molecule transitions in the lowest branch, which under¬ 
goes an impurity-trimer crossover before crossing with 
the molecular branch, there are two polaron to molecule 
transitions in the second-lowest branch. With increasing 
& c , the two transition points will move towards the BCS 
and the BEC limit, respectively. 

Similar picture for the polaron to molecule transition 
holds when we consider the special case where g ff = 
i.e., the fermion-fermion interaction and the fermion- 
impurity interaction are equal. As shown in Fig. [6^b), 
with an appropriately chosen cutoff momentum, the 
avoided crossing between the lowest two branches in the 
polaron spectrum is very narrow. If one starts from the 
ground state in the BCS limit and tune the interaction 
strength a s = = a& fast enough, the system would 

not end up in the trimer-like branch beyond the avoided 



- 

l 

bf-l 

Kf 

- 2 - 


(b) 


(k F a s ) x 0 


FIG. 6: (Color Online) (a) Polaron to molecule transition in 
the Q — 0 sector. The solid and the dashed lines correspond to 
the lowest two branches of the polaron energy, while the dash- 
dotted curve is the molecular energy relative to the molecular 
threshold. The fermion-fermion interaction is at resonance 
(as s= oo). Cutoff momentum is chosen to be k c = 10 k F , and 
the unit of energy is E F . (b) Polaron to molecule transition 
for the case of a s — = a^, with k c = 15 k F . 


crossing. In this case, there will be a polaron to molecule 
transition on the BEC side of the resonance, at which 
point the ground state of the system undergoes a first- 
order transition from the non-universal polaron-like state 
into a universal molecular state. This picture is consis¬ 
tent with the results reported in Ref. El- 


VIII. SUMMARY AND DISCUSSION 

We have studied the polaron excitations when an impu¬ 
rity is interacting with a Fermi superfluid. In particular, 
we show the importance of higher-order scattering pro¬ 
cesses caused by the superfluid fermions in evaluating the 
polaron energy. Consequently, the mean-field description 
of polarons becomes inadequate even for small fermion- 
impurity scattering lengths, especially when the fermions 
are tuned across resonance to the BEC side. This poses 
new challenges on the theoretical treatment of these sys¬ 
tems beyond the mean-field or perturbative approaches. 
Moreover, our work shows that the impurity-trimer cou¬ 
pling can be greatly enhanced by increasing the fermion- 
fermion and fermion-impurity interactions or by decreas¬ 
ing the effective range, which all lead to broader avoided 
level crossings. In practice, this can serve as a guideline 
to reduce three-body losses and maintain the stability of 
mixture systems. Our results can be directly probed in 
current cold atoms experiments. The polaron spectrum 
and its residue can be measured using radio-frequency 
spectroscopy mm®, and the stability of the system can 
be detected through the atom-loss measurement. 

Finally, we note that our polaron ansatz could be im¬ 
proved by including collective excitations in the Fermi 
superfluid due to phase fluctuations of the pairing field. 
This type of excitation dominates over pair breaking 
processes in the deep BEC regime of fermions, where 
the system can be modeled by a polaron on top of a 
molecular BEC [3QH38j , with impurity-boson and boson- 
boson scattering lengths respectively given by the three- 
body uni cm and the four-body J39I solutions. In 
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this work, we focus on the effects of pair breaking 
in the weak-coupling and the near-resonant regimes of 
fermions, where the collective excitations cannot quali¬ 
tatively change the main features of the polaron, i.e., the 
emergence of trimer physics, the avoided level crossings, 
and the non-universality of polaron spectra. 


manuscript. This work is supported by NFRP 
(2011CB921200, 2011CBA00200), NNSF (60921091), 
NSFC (11374177,11374283). XC acknowledges support 
from programs of Chinese Academy of Sciences. WY 
acknowledges support from the “Strategic Priority 
Research Program(B)” of the Chinese Academy of 
Sciences, Grant No. XDB01030200. 


Acknowledgements 

We thank Tin-Lun Ho, Hui Zhai, Yusuke Nishida 
and Xiwen Guan for helpful feedback on the 


[1] A. Schirotzek, C.-H. Wu, A. Sommer, and M. W. Zwier- 
lein, Phys. Rev. Lett. 102 , 230402 (2009). 

[2] S. Nascimbene, N. Navon, K. J. Jiang, L. Tarruell, M. 
Teichmann, J. McKeever, F. Chevy, and C. Salomon, 
Phys. Rev. Lett. 103, 170402 (2009). 

[3] C. Kohstall, M. Zaccanti, M. Jag, A. Trenkwalder, P. 
Massignan, G. M. Bruun, F. Schreck, R. Grimm, Nature 
485, 615 (2012). 

[4] M. Koschorreck, D. Pertot, E. Vogt, B. Frolich, M. Feld, 
M. Kohl, Nature 485, 619 (2012). 

[5] F. Chevy, Phys. Rev. A 74, 063628 (2006). 

[6] C. Lobo, A. Recati, S. Giorgini, and S. Stringari, Phys. 
Rev. Lett. 97, 200403 (2006). 

[7] R. Combescot, A. Recati, C. Lobo, and F. Chevy, Phys. 
Rev. Lett. 98, 180402 (2007). 

[8] R. Combescot and S. Giraud, Phys. Rev. Lett. 101 , 
050404 (2008). 

[9] N. V. Prokofev and B. V. Svistunov, Phys. Rev. B 77, 
125101 (2008). 

[10] M. Punk, P. T. Dumitrescu, and W. Zwerger, Phys. Rev. 
A 80, 053605 (2009). 

[11] X. Cui and H. Zhai, Phys. Rev. A 81, 041602(R) (2010). 

[12] S. Pilati, G. Bertaina, S. Giorgini, and M. Troyer, Phys. 
Rev. Lett. 105, 030405 (2010). 

[13] P. Massignan and G. M. Bruun, Eur. Phys. J. D 65, 83 

( 2011 ). 

[14] Y. Nishida, Phys. Rev. Lett. 114, 115302 (2015). 

[15] I. Ferrier-Barbut, M. Delehaye, S. Laurent, A. T. Grier, 
M. Pierce, B. S. Rem, F. Chevy, and C. Salomon, Science 
345, 1035 (2014). 

[16] X. Cui, Phys. Rev. A 90, 041603(R) (2014). 

[17] R. Zhang, W. Zhang, H. Zhai, P. Zhang, Phys. Rev. A 
90, 063614 (2014). 

[18] C. J. Pethick and H. Smith, Bose-Einstein condensation 
in dilute gases , Cambridge University Press (2nd edition, 
2008). 

[19] M. P. Marder, Condensed Matter Physics , John Wiley V 
Sons (2nd edition, 2015). 

[20] F. Chevy and C. Mora, Rep. Prog. Phys. 73, 11 (2010). 

[21] P. Massignan, M. Zcaccanti, and G. M. Bruun, Rep. 
Prog. Phys. 77, 034401 (2014). 

[22] While this is adequate to provide quantitatively accu¬ 
rate results in polarized Fermi systems when compared 


with results from quantum Monte Carlo studies EHE], 
we leave the justification of such a practice in the current 
system to future studies. 

[23] Y. Nishida, Phys. Rev. A 79, 013629 (2009). 

[24] D. J. MacNeill and F. Zhou, Phys. Rev. Lett. 106, 145301 

( 2011 ). 

[25] C. J. M. Mathy, M. M. Parish, and D. A. Huse, Phys. 
Rev. Lett. 106, 166404 (2011). 

[26] N. G. Nygaard and N. T. Zinner, New J. Phys. 16, 023026 
(2014). 

[27] S. Endo and M. Ueda, arXiv: 1309.7797 

[28] The residue here refers to the weight of bare impurity 
term in the background of the ground state Fermi su¬ 
perfluid. Thus, it does not include the excitation term 
with k + k x = 0, though the impurity still stays at initial 
momentum Q. 

[29] Close to (kFd r) - 1 ~ —4, an avoided crossing exists be¬ 
tween the second lowest and the third lowest branch. This 
is why the residue Z in the second lowest branch does not 
immediately go to small value (~ 0) beyond the avoided 
crossing with the lowest branch (see Fig. [TJbl)). 

[30] G. E. Astrakharchik and L. P. Pitaevskii, Phys. Rev. A 
70, 013608 (2004). 

[31] R. M. Kalas and D. Blume, Phys. Rev. A 73, 043608 
(2006). 

[32] F. M. Cucchietti and E. Timmermans, Phys. Rev. Lett. 
96, 210401 (2006). 

[33] M. Bruderer, W. Bao, and D. Jaksch, Europhys. Lett. 
82, 30004 (2008). 

[34] J. Tempere, W. Casteels, M. K. Oberthaler, S. Knoop, 
E. Timmermans, and J. T. Devreese, Phys. Rev. B 80, 
184504 (2009). 

[35] S. P. Rath and R. Schmidt, Phys. Rev. A 88, 053632 
(2013). 

[36] W. Li and S. Das Sarma, Phys. Rev. A 90, 013618 (2014). 

[37] F. Grusdt, Y. E. Shchadilova, A. N. Rubtsov and E. Dern¬ 
ier, arXiv: 1410.2203; Y. E. Shchadilova, F. Grusdt, A. N. 
Rubtsov, and E. Dernier, arXiv:1410.5691 

[38] R. S. Christensen, J. Levinsen, G. M. Bruun, arXiv: 
1503.06979, 

[39] D. S. Petrov, C. Salomon and G.V. Shlyapnikov, Phys. 
Rev. Lett. 93, 090404 (2004). 



